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ABSTRACT 

We have performed an extensive search for planet candidates in the publicly available 
K epler Long Cadence data from quarters Ql through Q6. The search method consists 
of initial de-trending of the data, applying the trend filtering algorithm, searching for 
transit signals with the Box Least Squares fitting method in three frequency domains, 
visual inspection of the potential transit candidates, and in-depth analysis of the 
shortlisted candidates. In this paper we present 150 new periodic planet candidates 
and 7 single transit events, 72 of which are in multiple systems. The periods of these 
planet candidates vary from ~0.17 day to ^440 da y. 124 of the planet ca ndidates have 
radii smaller than 3i?®. We recover 82.5% of the Batalha et al. (2012) KOI catalog. 
We also report 40 newly identified false positives- — systems that look like transiting 
planets, but are probably due to blended eclipsing binaries. Our search improves the 
statistics in the short period and small planet radii parameter ranges. 
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1 INTRODUCTION 

The field of transiting extrasolar planets (TEPs) has ex- 
ploded over the past few years. One major contributor is the 
Kepler Space Missi on, continuously mon i toring 156000 star s 



in a 115 deg field (Borucki et al 



2010 



Koch et al.||2010[ ). 



Since its launch in 2009, it has found 2321 planetary tran - 
sit candidate j 1 ] ( |Borucki et al.||2011a Batalha et al. 2012 1 . 
The first 4 months of data (Ql and Q2) yielded 1235 planet 
candidates associated with 997 host stars, including 60 con- 
firmed planets around 33 stars ( |Borucki et al.||2011b| here- 
after Bll). The Kepler team has developed increasingly so- 



phisticated procedures to identify planet candidates ( Smith 



et al. 112012 1 and multiple systems (Steffen et al. 112012 Ford 



et al. | 2012 1. Using more data (16 months) and improved 



detection procedures (the Transit Planet Search (TPS) al- 
gorithm; Jenkins et al.|[2010c Tenenbaum et aT7||2012 1, the 
total number of planet candidates has almost doubled since 
the release by Bll ( jBatalha etliL|2012[ hereafter B12). 

An independent search with different tools can build 
confidence in the reliability of the detections, and the com- 
pleteness of the sample, by e.g. providing new candidates 
that were missed by the Kepler science team. One such ef- 
fort is the citizen science initiative, called PlanetHunters 

based on the idea 



(Fischer et al. 2012 Lintott et al 



2012) 



of Zooniverse ( |Lintott et al.||2008 l. Making use of human 



* E-mail:xuhuang@princeton.cdu; 

1 http:/ /kepler. nasa.gov/Mission/discoveries/candidatcs/ 



eyes to search for transit-like events through a user friendly 
computer interface, 6 new planet candidates were published 
by far ( jFischer et "aX|2012| |Lintott et aL]|2012[ ). These were 
then subjected to the vetting procedure of the Kepler team, 
and five of them survived this process, i.e. they are not false 
alarms, as much as it can be determined from the Kepler 
data. However, conducting a visual search of all the raw light 
curves takes a lot of human effort. Quoted from |Fischer et| 
al. (2012 I: "It is impractical for a single individual to review 
150, 000 light curves in every quarterly release 
While challenging, it is, however, 
150, 000 light curves, 



each of the 

of the Kepler data base.' 
feasible for an individual to examine ■ 
if sophisticated computer algorithms narrow down the list 
to a somewhat smaller sample of candidate transits that can 
be then checked very carefully. Also, note that small plan- 
ets are often hard to recover by visual inspection without 
phase-folding the light curve at the suspected periodicity of 
the signal. 

Encouraged by these observations, we started an inde- 
pendent search for transiting planet candidates in Kepler's 
long-cadence data. Kepler observations are grouped into so- 
called quarters (each having a duration of 3 months, ex- 
cept for Q0 and Ql, which are shorter). At the end of each 
quarter, the spacecraft rotates 90 degrees, to adjust its solar 
panels. The majority of Kepler stars are observed as long 
cadence (LC) targets, for which data are obtained by gath- 
ering over 270 exposures within 29.4 min; 512 stars are also 
observed in short cadence (SC) mode with 58.9 s intervals 
( |Jenkins et "aT1|2010b| |Gilliland et~aT1|2010| |Murphy|[2012 1 . 



© 0000 RAS 



2 X.Huang, G.A.Bakos, J.D.Hartman 



We used quarters Q0-Q6 in this search, which data were 
released to the public in Jan 2012. 

Our methodology is based on our experience conduct- 
ing a similar search with HATNet (Bakos et al. 20041, a 



wide-field ground-based survey. Broadly speaking, the Ke- 
pler space-based data is of much higher quality than any 
ground-based data. Ground based observations often exhibit 
large gaps in the time-series, either due to the rotation of 
the Earth, or inclement weather conditions. The data qual- 
ity is highly variable due to changing extinction, clouds, 
background, seeing, etc. The per-point photometric preci- 
sion is typically worse than from space, partly because of the 
above effects of ground-based observations, and partly due 
to the use of inexpensive hardware (e.g. front illuminated 
CCDs). Altogether, our ground-based data is of lower signal- 
to-noise, has inhomogeneous quality, and exhibits complex 
systematic variations with long gaps. For example, not a 
single transit event has been found in HATNet data by di- 
rect visual inspection; transits are detected through sophis- 
ticated data mining and phase-folding. Also, no robust sin- 
gle transit event was ever found by HATNet. Consequently, 
tools developed for a transit search using ground based data 
may be very efficient in recovering transit signals from Ke- 
pler. 

In this paper we employ tools from the transit detec- 
tion pipeline of the HATNet project, after sufficient modi- 
fications to conform to the Kepler data. We also develop a 
pre-filtering method that corrects the Kepler light curves for 
known anomalies and systematics, before searching them for 
transit events. Our search is blind in the sense that the list 
of Kepler candidates was not consulted during the search. 

The structure of the paper is constructed as follows. The 
data processing methodology is discussed in Sj2] We laid out 
our findings in ^3] including the recovery of Kepler planet 
catalog, our selection and modeling of new candidates, and 
the properties of these new candidates. We make our con- 
cluding remarks in Q 



2 DATA ANALYSIS 

2.1 Removal of points and long trend filtering 

We make use of the Kepler public LC light curves. We begin 
with the TIME and SAP_FLUX (raw flux) columns in the 
FITS files. The raw flux is already corrected for cosmic rays 
and background variations by the Kepler team ( |Jenkins et| 
al.|2010af First we convert the fluxes to magnitudes and set 
the mean value for each star to its Kepler magnitude taken 
from the Kepler Input Catalog (KICQ 

The second step is to clean the light curves based on 
the data anomalies table in the Kepler Data Release Notes 
for each quarter ( |Christiansen et al.||2012|b| |Machalek et| 



al. 2010 20111. We summarize all the important events in 



Ql through Q6 in Table [T] We also describe the definition 
of anomaly types (adopted from Kepler Data Characteristic 
Hand Book) and our methods of correction for each type in 
the table notations. Generally, for data anomalies involving 
a discontinuity we model the jump by a polynomial with an 
offset in magnitude after the gap. The fitting uses 50 points 



http:/ /tdc- www. harvard.edu/softwarc/catalogs/kic. html 



on both sides around the gap. The offset is then subtracted 
from the data after the jump. In the time range during safe 
modes (Q2 and Q4), or earth point recoveries (Q3, Q4 and 
Q6), there is an exponential decay in the flux. We iden- 
tify and remove the whole exponential decay rather than 
attempting to correct for the effect. 

Most of the stars exhibit long-term trends. We follow 
the traditional high-pass filter method ( Ahmed et al. |1974 



Mazeh & Faiglcr 20101, and apply a cosine filter on all the 



cleaned light curves. For each light curve, before computing 
the filter, we generate a model by applying a 100 point (~ 2 
day long) median filter. This is done to prevent distortions 
due to outlier points (including introducing spurious "tran- 
sit" signals). The cosine filter is then computed as the sum 
of a linear component and N = T tota i/AT cosine functions, 
where the highest frequency is 1/AT (AT = lday), and 
Ttotai is the total time span of the light curve: 



M(t 3 ) = a 



(tj-ts) 
Ttotai 



+ 



E 



Ttotai 



(i) 



Here tj is the time of the j'th measurement, and t s is the 
first time instance in the light curve. The coefficient a for 
the linear component and coefficient bi for the ith cosine 
functions are computed by a least squares fitting procedure 
on the model. The fitted trend M(tj) is then subtracted from 
the light curves. We apply this cosine filter to light curves in 
every Quarter separately, and then combine the long trend 
filtered light curves from Q1-Q6 by offsetting the magnitude 
of all the quarters referring to the magnitude of Ql. 



2.2 Sky groups and Systematic Trend filtering 

Following the long-term trend filtering procedure described 
above (§ |2.1[ ), we then apply the Trend Filtering Algorithm 
(TFA) developed by Kovacs et al. ( 2005 1 on the combined 
(Q1-Q6) light curves. The idea of TFA is to select a set of 
template light curves, which we assume to contain informa- 
tion of the systematic variations, and then construct a linear 
filter based on their shared time-series for each light curve to 
be corrected. TFA can remove systematic variations that are 
either shorter time-scale than those corrected by the cosine 
filter, or have an arbitrary functional form that is not well 
described by the sum of cosine functions. TFA assumes that 
the light curves are sampled at the same time-instances. 

To construct a set of template light curves with the 
same time-base, we make use of the sky group information 
provided by the Kepler team. The sky group number is de- 
fined as the CCD channel number on which the stars fell 
during Q2 of the Kepler operation. Kepler has 21 modules. 
Each module contains 2 CCD chips and each CCD contains 
2 output channels. The focal plain rotates 90 degrees when 
the telescope switches to a new quarter every 3 months, 
except for the initial transition between quarters Q0 to Ql. 
Generally, the stars that belong to the same sky group share 
the same CCD channel all the time, although this channel 
changes from time to time. Therefore, each sky group shares 
the same time base, and has similar instrumental systematic 
trends. Additionally, stars in the same sky group are related 
in terms of their sky position, so that the local variabilities 
(e.g. local background variations or flux contamination from 
nearby stars) are often shared. TFA is designed to reduce 
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this kind of general (shared by a number of stars) system- 
atics. 

We construct a separate filter for each sky group, se- 
lecting ~ 300 template light curves in each case. Template 
stars are selected randomly, but in a manner that ensures 
a uniform distribution of their positions across the field. To 
exclude variable stars, we also impose a constraint on their 
median deviation around the median magnitudes (MAD, 
a quantity that is insensitive to outliers); stars with high 
MAD are ruled out from the templates. The total number 
of data points (time samples) in each template time base is 
~ 17000 — 22000, depending on the total length of time se- 
ries in the sky group. In other words, we are not overfitting 
the light curves. 

Altogether 124840 stars in 84 sky groups were selected 
and analyzed with TFA from the Kepler public LC data. 
For each sky group, we only selected stars that have been 
observed during the complete time range, and have not been 
affected by the failure of module 3. This is because the TFA 
analysis requires the same time-base to generate the filter 
per sky group. Module 3 failed in the middle of Q4, while 
observing sky groups 5, 6, 7 and 8. Due to the rotation of 
spacecraft, the Q5 data in sky groups 49, 50, 51 and 52, and 
the Q6 data in sky groups 77, 78, 79, and 80 were not avail- 
able. Stars in these sky groups are still included if they have 
a complete data set in other quarters. Otherwise, we only 
selected light curves containing all the Q1-Q6 data (which 
has a total observation time of ~ 500 days). 

We note that the Kepler team has recently implemented 
a new "cotrending" algorithm, called PDC-MAP ( |Smith et] 
aT] |2012[ ). This algorithm uses 16 Cotrending Basis Vectors 
(CBV) that arc generated by a singular value decomposition 
method applied separately to each channel and each quarter. 
This method has some commonalities with TFA. We did not 
make use of the PDC-MAP data, since it was not available 
for the Q1-Q6 data when we started our analysis. 

In Figure [T] we demonstrate our filtering process on a 
randomly selected light curve (KIC 003346154). The great 
improvement by the procedure is clearly demonstrated. 

2.3 Box Least Square Fitting (BLS) and Transit 
Analysis 



We use the box least square fitting algorithm (Kovacs et 



al. 2002 \ to search for periodic transit signals in the TFA- 



filtered light curves. In order to maximize the efficiency of 
the BLS search in a wide frequency range, we divide the fre- 
quency search into three, only slightly overlapping frequency 
domains: 0.3 d- 1 < /i < 9.0 d _1 , 0.02 d" 1 < f 2 < 1.0 d" 1 
and 0.005 d" 1 < f 3 < 0.03 d" 1 . We use a different number 
of frequency steps and BLS bins in each domain. We use 
the SNR (signal to noise ratio, measured in the BLS spec- 
trum) and DSP (dip significance) parameters (for details, 
Kovacs et al. 2002) of the first five frequency peaks re- 



ported by BLS to select candidate transit signals for manual 
inspection. 

We adopt selection threshold of 11 and 8.5 for SNR 
and DSP in the middle and long period range. The selection 
threshold of short period range is somewhat higher (30 and 
20 for SNR and DSP) considering that BLS tends to respond 
to short period easier and there are fewer KOIs in this pe- 
riod range as our reference. On top of these requirements, 



we neglect the BLS peaks with frequency too close to the 
frequency domain boundary (i.e. the period ranges used in 
selection are 0.67 d" 1 < f[ < 6.67d _1 , 0.0202d" 1 < /_ < 
0.67d~ 1 and 0.005 d" 1 < f 3 < 0.0202 d" 1 ). We also reject 
those with a transit duration much longer than expected 
(qi < 0.5, (/2 < 0.1 and qg < 0.013) or a transit depth in- 
dicating very large planetary radius (dip > 0.04). These are 
not optimized selection criteria, but instead ensure that no 
shallow or rare transit events are missed. The low detection 
threshold also means that events will require close visual in- 
spection. With the above limits, we selected ~10% of the 
stars to fold with the BLS peaks that satisfy the above re- 
quirements. The number of BLS peaks selected is ~3% of all 
the first five best BLS peaks. We then manually inspected 
all the folded light curves. 

We take three random sky groups as an example to 
illustrate our selection process in Figure|2] We show the dis- 
tribution of BLS peaks (the mid period range) in the SNR- 
DSP plane. The solid lines are our lower limits for selection. 
The green dots are selected periods for further examination. 
The best BLS peak of KOI planet candidates are represented 
with red crosses. Majority of the KOIs are selected with high 
SNR and DSP. 

Transit-like features in the folded light curves are 
flagged and further examined in the visual inspection. We 
reject light curves with recognizable depth variation for odd 
and even peaks. We also check the harmonics of the de- 
tected period to ensure the detection of a correct period. 
For the transits visible in the unfolded light curves, we di- 
rectly inspect the shape and transit center of each transit; 
for the transits invisible in the unfolded light curves (due to 
low S/N), we perform phase-folding with the BLS-detected 
frequency before examining the data. 



3 RESULTS 

3.1 Comparison of Our Sample to that of the 
Kepler Team 

We flagged 2180 stars as possible transiting planet hosts 
during manual inspection. Among them, we found 18 stars 
which are categorized as eclipsing binaries (EB)FjdPrsa e ^ 
al]|201l| ) or false positives (FP) P] pBorucki et al.[2011b[ ) by 
the Kepler tea m. We also cross-matched our results with th e 
KOI catalog pi (Borucki et al.|2011a|b| |Batalha et al.|2012 |. 
We note that our detection efficiency is comparable to the 



TPS algorithm used by the Kepler team ( Tenenbaum et al 



|2012| Hereafter T12). Their algorithm yields 5392 detections 
and detected 88.1% of the 1235 KOIs in the BU catalog us- 
ing data from Q1-Q3. We recovered 92% of the BU catalog 
in our analysis (1124 of the KOIs in BU are in our initial 
data set, from which, 1034 are flagged) thanks to the longer 
time base we used compared to T12. 

We demonstrate the selection of KOIs compared to the 
B12 catalog in Figure [3] There are 1518 KOI stars (corre- 
sponding to 1982 KOI planet candidates) from the B12 cat- 
alog included in our overall LC light curve samples. The rest 



3 http://archive.stsci.edu/kepler/eclipsing_binaries.html 

4 http:/ /archive. stsci.edu/kepler/false_positives. html 

5 http:/ /archive. stsci.edu/kepler/planet_candidates. html 
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Figure 1. An example of our filtering process as applied to the light curve of KIC 003346154. The ordinate is in the units of approximate 
Kepler magnitude: Top panel: correction of safe mode (left) and tweak (right), (the example is taken from Q2). The corrected light curve 
is plotted below the uncorrected one. Middle panel: combined light curve from Q1-Q6 after the pre-filtering process, i.e. correction of 
anomalous data points and combining the light curves from all the quarters. Bottom panel: light curve after long trend filtering (solid 
symbols) and after TFA (empty circles). 



of the KOI stars in the B12 catalog either fail to fulfill our 



long cadence time baseline length requirement (see [ 2.2 1 or 
have transit depths greater than 0.04, and are therefore re- 
jected by our procedure. We found that 1311 (86.4%) of the 
KOI stars are flagged in our selection process. 1636 (82.5%) 
of the KOI planet candidates are detected (either with the 
correct period reported by BLS or are detected with the 
wrong period but later recovered during visual inspection). 

We failed to recover about 17.5% of the KOI planet can- 
didates, mostly due to no significant peaks found in the BLS 
spectra corresponding to the transit periods. This partly re- 
sults from the choice of BLS parameters and also the non- 



periodic properties of some transit signals. The hard cut 
in DSP is the second most important reason for the rejec- 
tion of some KOI planetary candidates. Some KOIs have low 
DSP either due to over correction of our light curves or light 
curves with high noise level. There are also 40 KOI planet 
candidates with BLS peaks selected by our pipeline but re- 
jected during visual inspection. From further investigation, 
we find that 32 of these rejected transit signals are not vis- 
ible when folded with twice the detected periods, probably 
due to strong transit timing variations. A detailed detection 
report of all the KOIs is presented in Table [2] We do not 
model the KOIs unless the host stars are found to have other 
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erated in a frequency range adjusted according to the visual 
determination of the frequency of the stellar variation. To 
preserve the transit signal, we applied a median filter with a 
window width at least twice that of the detected transit du- 
ration, before generating the cosine filter. BLS was applied 
again on the corrected light curves. The first ten peaks in 
the BLS spectrum are examined and compared to the previ- 
ous detection, ensuring the detection of the transit signal is 
robust. As shown below in j ]3.3[ all candidates were checked 
against false positive scenarios. 



10 100 
SNR 

Figure 2. A demonstration of our selection methods. SNR versus 
DSP for all the first five BLS peaks from stars in sky groups 1, 
12, and 19 (black dots). Solid lines present the lowest limit of 
SNR and DSP for selection. Green asteroids resemble all the BLS 
peaks actually selected for folding considering the period range 
and transit duration. Red crosses show the most significant peak 
of all the KOI planet candidates (34 in total) in these sky groups. 
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Figure 3. Summary of our recovery of KOIs in the B12 catalog. 
We show the number of stars (S) and planet candidates (P) in 
each category. The recovery rate of KOI host stars in our search 
range is 86.4%. The recovery rate of KOI planet candidates is 
82.5%. For more information, see jpTTland Table [2] 

high SNR and DSP transit signals that are not reported by 
the Kepler team. The additional signals are treated as new 
candidates in our short list. 

3.2 Generation of the planet candidate short list 

We examined in detail the remaining candidates which are 
not included in the publicly available lists of candidates, 
EBs, or false positives (FPs). We re-applied the cosine filter 
to these light curves constructed with a greater number of 
cosine functions and a smaller frequency interval. We aimed 
to clean systematic variations from the light curves (whether 
intrinsic to the stars, or due to instrumental effects), leaving 
flat light curves with transits. The cosine functions are gen- 



3.3 False positive detection and robustness 
checking 

We use the moment-derived centroids provided by the Ke- 
pler FITS files to further eliminate possible FPs. For long 
period planet candidates, we visually examine the centroids 
at the transit time. For short period planet candidates, we 
examine the phase-folded centroid curves (de-trended by the 
cosine filter). We present one of our "failed" transit can- 
didates as an example in Figure [4] The folded light curve 
would naively suggest that there is a transit signal due 
to a planet with Hp/Ft* = 0.0198 ± 0.0017 and period of 
P = 2.00783T3.8 x 10~ 5 day. However, the phased centroids 
show a shift of ~ 0.004 pixels in the y direction during tran- 
sit events, which indicates this is likely to be a false positive 
signal due to a blended eclipsing binary. According to the 
2MASS image stamp, there is no nearby source within an 
area of 20" x 20", i.e. the binary is not resolved in the 2MASS 
images. We provide a list of forty false positives flagged by 
this method in Table [3] These stars have not been reported 
by the Kepler team in their false positive lists. They are also 
not reported as candidates by the Kepler team. The stellar 
information and estimated shifts in both directions are also 
given in our table. 

We also used the public target pixel files and PyKE 
packag^] developed by the Kepler team, to obtain pixel light 
curves for all of our candidates with transit depths greater 
than 1 mmag. It is difficult to perform the same analysis on 
shallower transits, because of the low signal-to-noise of the 
events in the light curves of individual pixels. 

Here we take one of our new detections, KIC 005437945, 
as an example for our photometry analysis. We detected 4 
non-periodic transit events altogether, which could be ex- 
plained as due to two long period planet candidates with 
different epochs, depths and durations. The transits in Ql 
and Q6 are due to a P w 440 days planet candidate; and 
the transits in Q2 and Q5 are due to a P « 220 days planet 
candidate (i.e. they are in a 2:1 resonance). We present the 
analysis for the transit event in Q6 here. We show the pixel 
image of out-of-transit, in-transit and the difference imag- 
ing in Figure [5] The images are computed by plotting the 
mean out-of-transit flux (±2 day around the transit events), 
the mean in-transit flux and the difference between the two. 
The difference imaging during transit is identical to the out- 
of-transit flux distribution. No obvious background source is 
indicated for this particular star. 

We further use the pixel calibration technique enabled 
by PyKE to extract light curves from every single pixel in 

6 http:/ /keplergo. arc. nasa.gov/ContributedSoftwarePyKEP.shtml 
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Figure 5. Out-of-transit image (in log scale of flux), in-transit image (in log scale of flux) and the difference between the two (in linear 
scale of flux) from the pixel files for K1C 005437945 Q6 transit. We do not see a visible shift in the flux distribution on the pixel image 
during transit, demonstrating that the apparent transit is not due to a variation of a background source. 
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Figure 4. Light curves and centroids for KIC 004072333, phase- 
folded using the detected period and epoch. Top panel: phase- 
folded light curve (bottom) with period P = 2.00783 days and 
E = 2455302.8774 (BJD) and its residual from the fitting (top). 
Bottom panel: folded centroids in the x direction (top and grey) 
and the y direction (bottom and dark) with the same parameters. 
This detection is flagged as a false positive due to the shift of 
centroids in the y direction during transit. There is no visible 
companion in the 2MASS image stamp within 20" X 20" for this 
particular star. 



the aperture. Figure [5] shows the light curves extracted sep- 
arately from 6 pixels in the Q6 aperture (the aperture is 
shown in Figure We can see that while the magnitude in 
every single pixel changes during the transit, they all show 
visible evidence for the transit event with roughly the same 
depth. The averaged flux from all the pixels has a flat out-of- 
transit magnitude. In addition, the centroids do not present 
an anomalous shift during the transit events. 



For the rest of the transits which are not suitable for ap- 
plying the single pixel analysis, we use other methods to at 
least ensure the robustness of the detection. We removed the 
detected transit signal according to the epoch, period and 
transit duration reported by BLS. Then we computed six 
different BLS spectra for each of the new light curves. The 
BLS spectra were obtained with 3 different sets of parame- 
ters (different frequency steps and bin numbers), each set of 
parameters are used twice. The first 10 peaks of every BLS 
spectrum were analyzed. We compared these 60 frequency 
peaks to our original transit detections. For all of the plane- 
tary transit signals we analyzed, the original periods as well 
as their harmonics were not detected in the transit-removed 
light curves. 



Altogether 150 planet candidates entered our final list 
of new detections. The properties of their host stars are pro- 
vided in Table[4] The detected periods range from ~ 0. 17 day 
to ~ 440 days. We also found 6 new multiple systems (alto- 
gether with 15 planet candidates) in stars not in the B12 
catalog. We found 57 new transit signals in KOI hosting 
stars; 43 of these are also independently reported by |Ofir| 
& Dreizler (2012 Hereafter 012). We also include 7 single 
transit events. One of these single transit events is around 
a KOI star with a known planet candidate. Using the con- 



vention of Batalha et al. (20121 for single transit events, we 



assign a negative integer period number for these potential 
candidates. We compute the minimum allowed period ac- 
cording to the given time span of the light curve and the 
epoch of transit. This is taken as the estimated period for 
single transit events in modeling. 
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Figure 6. Red solid line: light curve flux for transit events in 
Q6 for KIC 005437945. Blue triangles: raw centroids x (top panel) 
and y (bottom panel). Different types of black points: light curves 
generated from different pixels separately in the aperture using 
kepextract from the PyKE package. 



3.4 Analysis 

The transit modeling of all the candidates is based only on 
the Kepler light curves and the stellar parameters provided 
by the KIC. The parameters we can obtain directly from 
light curves are the transit depths, durations, ingress/egress 
durations, and individual transit centers. 

Without radial velocity (RV) data, mass determinations 
for these systems are generally not available. It may be 
possible to measure the masses through subtle photomet- 
ric effects, like ellipsoidal variations and relativistic beam- 
ing ( |Mazeh fc Faigler|2010||Kipping fe Spiegel|20lTj ). These 
effects are prominent for close binary stellar systems as well 
as massive planet companions (Loeb & Gaudi 2003]) . We did 
not observe these effects in our candidates. 

The eccentricity is also unknown, although broad lim- 
its can be placed on the orbital configuration for very wide 
transits (Kipping 2008). An eccentric orbit could result in 
asymmetry in the transit light curves, as well as a shift in 
the mid-time of transits relative to the occultation events. In 
principle, the eccentricity could also be derived from model- 
ing the detailed shape of a transit light curve. However, de- 
tecting these effects requires extremely high resolution and 
SNR light curves, which were not available for our candi- 
dates. Generally speaking, we observe no apparent asym- 
metry in any of our candidates, which suggests modest ec- 
centricity or certain values of argument of periastron. We 
assume circular orbits in our modeling, following the con- 



vention of previous KOI modeling with only transit light 
curve information ( jBatalha et al.|2012 l. 

We assume no flux dilution from blended nearby stars 
in our modeling. This is plausible for most of our candi- 
dates. We inspect the 2MASS image stamps with an area of 
20" x 20" centered on each target. Four of them have nearby 
companions, which are marked out in Table [5] 

In theory, one could also fit for the limb darkening coef- 
ficients, but this requires very high signal-to-noise and well 



sampled data, i.e. deep transits or many transit events ( Kip 
|ping fc Bakos 2011| ). We use a quadratic limb darkening for 
malism to model the transit light curves. The limb darkening 
coefficients (LDC) corresponding to the stellar atmospheres 
are interpolated with the stellar parameters from KIC in 
the ATLAS model grid LDCs provided by |Sing| ( |2010[ ) for 
Kepler. The stellar parameters and the LDCs for the candi- 
dates are listed in Table [4] We also generate a grid of 1000 
randomly selected stars from KIC with known stellar param- 
eters and the interpolated their LDCs following the method 
described above. We then use these to determine LDCs for 
stars without information such as T e s, logy and [Fe/H] by 
linear interpolating in the J, H, K color space. Since the 
LDCs obtained with this method have larger uncertainties, 
the derived planet parameters for planets without host star 
parameters provided by KIC should be treated with caution. 
To conclude, in our transit modeling, we do not fit for: 

a) the planet mass M p ; 

b) the eccentricity e of the orbit, which is assumed to be 
zero; 

c) the blending due to nearby stars, which is set to be zero; 

d) the LDCs of stars, which are computed as described 
above. 

We did fit the following geometric parameters which 
correspond: 

a) the fractional planet radius R p /R,; 

b) the square of the impact parameter square b 2 ; 

c) the inverse of half duration (/R, — 2/Tdur- This quantity 
is related to a/ R* for zero eccentricity via the relation: 

a 2tt 1 



R* 



R, p ^/T^W 



(2) 



In the fitting procedure we constrained the quantities 
b 2 ,R r / R*,C/ R* to be the same for each individual transit 
for a selected candidate. 

We also fitted for the additional parameters of out-of- 
transit magnitudes and transit centers. For the short period 
cases (with a period shorter than 30 days), we fit the median 
out-of-transit magnitude, first transit center Ta and the last 
transit center Tb as free parameters, with the total number 
of transits N fixed. We also assume that the transits are 
strictly periodic. We do not model the out-of-transit varia- 
tion; the light curve is assumed to be flat. For our long pe- 
riod candidates, we use a slightly different method. We take 
the out-of-transit magnitude and transit centers of individ- 
ual transits as independent parameters, the total number of 
free parameters in the fit for a transit light curve with N 
transits is 2iV-r-3. The periods of the candidates with only a 
single transit were estimated by the stellar density through 
KIC parameters, and by assuming no limb darkening and 
circular orbits (Yee & Gaudi 2008). The estimated periods 



for these cases are shorter than the lower limit constrained 
by the fact that only one transit event was observed during 
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Q1-Q6. These single events were fitted as if they were an in- 
dividual transit in a long period system with a rough lower 
limit on the period. 

In the fitting we used the formalism of Mandcl & Agol 
( 2002 1 , and the methodology laid out in the analysis of H AT- 
Net planet discoveries (Bakos et al. 20101. A 10000 step 
Markov Chain Monte Carlo (MCMC) simulation is then ap- 
plied around the best fitting parameters to explore the pa- 
rameter space. The final reported planet parameters and 
estimated errors are taken as the median and median devi- 
ation of all the accepted jumps in the chain. The period is 
then recalculated by taking the median of {Tb — Ta)/N for 
all the accepted jumps in the chain. We derive the transit 
number closest to the average of Ta and Tb (weighted by 
their errors as derived from the MCMC runs), and use the 
transit center of this event (calculated from Ta, Tb, N) as 
the optimal epoch. 
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Figure 7. New Long Period planet candidates. For each candidate we show the best-fit model on the phase-folded light curve (bottom 
of the panel) and the residuals (top of the panel). For each candidate, we model all transits separately. We then fold the transits using 
the locally computed transit centers. We show the KIC number, period and scales of the figure of every planet candidate at the top of 
the subfigures. The x, y scales of the subfigure size are marked on the top right, in the units of [hours, lmmag] . 
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Figure 8. (a) 



© 0000 RAS, MNRAS 000, 000-000 



150 New transiting planet candidates from Kepler Q1-Q6 data 




— i 1 1 1 1 — - 

K I ™»746 1 3. 1 '= I I >hv) 1 1 2,.n.: 



—\ — p — i — i — i — i — i — 

. KKflfc J [ 953 b, ]fe [ .'J9(]iny] | S.4.5J 



i i i i i i i 

KIC...f4:^Mii.l 1 =.^.i!viUiyi |lr..4| 



'tis > • ■>£%' 



t 1 — i — n — I — r — i — n — 

- KKEQ5542SSJ, e=4..77(Day) [a,l.fl] 




3 VI 




- KlftXWfiSSfiSS, P=2.79(Day) 16,21 



Kl( DD57 7} Ul , EMLSSfDay) , | p,3 JL] 



_L I ,-L. 



_L, 1_ 




K\axw'n$\. i 3 =K..^(HEiy) [12,4] 




H 



Kit "Oa^SS I 9. I 3 =S. 24( I lay) 1 20.4.5 1- - 



~i "n 1 r 1 — 

KIO1060+7072, P=4.27(Day) [lQ^.S] 



is ; 



H 1 1 h 



Kl £ (KIM.') W.K P=.l. fi'Ji I Itiy i I S.-l 1 







M c (hour) 



t-t (hour) 



t-t c (hour) 



Figure 8. (b) 
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Figure 8. (c) 
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Figure 8. (d) 
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Figure 8. Short period planet candidates. For each candidate we show the best-fit model on the phase-folded light curve. The residuals 
are flat so we don't show them in the figures. The x, y scales of the sub-figure size are marked on the top right, in the units of [hours, lmmag]. 
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Figure 9. Single transit events. For each candidate we show the best-fit model on the light curve, together with the residuals from 
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[hours, lmmag]. 
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3.5 Properties of candidates 

We summarize the planet parameters for all the candidates 
in Table [5] The best fitted model on the phase folded light 
curves and the residuals from the best fitted model are pre- 
sented in Figures [7] [8] and [9] for long period, short period 
and single transit planet candidates. 

We compare the planet candidates reported in this pa- 
per with the B12 catalog KOIs in Figure [l0| Our new planet 
candidates follow a similar distribution in the SNR and DSP 
space as the KOIs. This demonstrates that our recovery of 
new candidates is due to the efficiency of our pipeline rather 
than due to a reduction in the selection thresholds. We also 
show in the right of Figure [10] that the majority of new 
candidates have short period (less than 10 days) and small 
transit depth. For comparison, we plot the KOIs missed by 
our search in green, suggesting that our detection sensitivity 
drops as the period of the transit signal increases. 

We plot the period and planet radius (in units of host 
star radii) histogram in Figure [TTj The turn over points in 
both the period and radius space are modified with the sup- 
plementation of our new candidates. We suggest that our 
methods are more sensitive to short period candidates since 
the detection efficiency of BLS is generally higher in high fre- 
quency. We also confirm that the KOI samples are almost 
complete in the mid-period range. We do not have much 
advantage in the long period range over the original Kepler 
pipeline. The sample size around P ~ 100 days slightly gains 
with our new candidates. Signals with periods longer than 
100 days usually suffer from over-correction of individual 
transit or low DSP simply because of lack of points in the 
transit phase. 

Finally, we provide remarks on some of the interesting 
systems below Q 

KIC 005185897: We found three planet candidates in 
the system, all with transit depths ~ O.lmmag. The periods 
of the candidates are 4.08 days (a), 11.42 days (b) and 6.83 
days (c). (See Figure [8] (a) , row 6, column 2, 3 and 4.) The 
stellar radius from KIC is only 0.51 Rq, which makes the 
the modelled radii extremely small for all the candidates in 
the system. The P = 11.42 day signal and 4.08 day signal 
are modelled to be ~0.62 R®, and the 6.83 days signal is 
slightly larger (0.67-Re). They are also the smallest planet 
candidates we found around non-KOI stars. Both the pair b 
and c, c and a are around 5:3 resonance. The planet candi- 
date (a) was also found by the TPS algorithm as a potential 
transit signal (T12). 

KIC 005437945: We identify four transit events in the 
light curve. The transit events in Ql and Q6, and the ones 
in Q2 and Q5 share the same depths and durations, respec- 
tively. The former is modelled to be a 7.5 R® planet can- 
didates with 440 day period; the latter is modelled to be 
a 6.4 T?^ planet candidate with a period of 220 days. (See 
Figure [7| row 2, column 2 and 3.) We note that the transit 
in Ql is independently identified and classified as a single 
transit event by the Planet Hunters ( Schwamb et al.|2012 |. 
A tweak happened during the transit in Q2 for the inner 
planet candidate; by carefully offsetting the magnitude on 




Figure 12. Image of the star KIC 5437945 taken with the APO 
3.5 m Echelle slitviewer. The guider is centered on the star. A 
companion is resolved within ~ 5". 
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Figure 13. The two planet candidates of star KIC 5437954: 
from black to grey are the phased folded light curve, and the 
phased folded x and y direction centroids, all displayed separately. 
The upper (lower) three plots are for the first (second) planet 
candidate, with period P ~ 440(220) day. The centroids in the 
figure are the mean centroids after de-trending. 



The planet and stellar parameters we quote here are only ap- 
proximated numbers for easy comparison. For accurate parame- 
ters and estimated errors, refer to Table [Bl 
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Figure 10. Comparison of KOIs and the planet candidates reported in this paper. Left: in SNR-DSP space; Right: in period-transit 
depth space. The black dots represents all the KOI planet candidates. The red cross represents the new candidates in this paper. In the 
right figure, the green squares show the KOIs missed by our searching process. 




both sides of the tweak, we recover the full transit signal. 
We resolved a faint close companion (~ 5") in the Echelle 



slitviewer of the APO 3.5 m telescope (See Figure 121. We 



examined 2MASS image stamps in J,H and K and DSS 1,2 
images in red and blue. All of the above show that the com- 
panion is bluer than KIC 005437945, which indicates that 
it is unlikely to be a physically associated companion. We 
present the phase folded de-trended centroid for both planet 
candidates in Figure (13] No anomalous motion is shown in 
either direction during transit. 

KIC 005965819: We found two planet candidates in 
the system, both with modelled radii less than 2 R® . The 
inner planet candidate orbits the host star with a period of 
8.2 days; the outer planet candidate has a 19.6 day period. 
(See Figure |8](b), row 3, column 1 and 2.) They are not in 
a tight resonance with each other. We note that the transit 
durations of these two planet candidates are longer than the 
values expected assuming circular orbits. This might indi- 
cate that the planets are in eccentric orbits or the stellar 
radius of KIC 005965819 is in fact larger than what is listed 
in the KIC. 

KIC 007673192: We found four planet candidates in 
the system, all with modelled radii around 1 . The periods 



of the candidates are 6.12 days (a), 4.03 days (b), 8.92 days 
(c) and 11.24 days (d). (See Figure [8] (c), row 3.) They are 
not around any low order resonance pairs. T12 also identified 
the P=6.12 day signal as a potential transit signal. 

KIC 009962455: We found two planet candidates in 
the system, both with modelled radii around 1 . The pe- 
riods of the candidates are 23.20 days (a), 5.78 days (b), in 
4:1 resonance. (See Figure [8] (d), row 5, column 2 and 3.) 
The 23.20 day component in the system was also found by 
T12. 

KIC 009535585: The shortest period candidate we 
found. The signal is modelled as a 1.6 R® super earth with 
a period of only 0.17 day. (See Figure [8jd), row 2, column 
2.) If this candidate is confirmed, it would hold the record 
of the shortest period among the KOIs. 

KIC 011253827: We found two planet candidates in 
the system. The outer planet is Saturn-like with a radius 
of 6.18 R® and period of ~ 88 days. The inner planet is 
~ 3.82 i?0 with a 43 days period. (See Figure [7] row 6, 
column 3 and 4.) The Planet Hunters independently found 
these signals in their web discussion. 

KIC 012692087: Also a extremely short period candi- 
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Figure 14. Folded de-trended light curve of KIC 005689351. 
Each transit is folded relative to its measured epoch, with the 
signal corresponding to the other planet candidates removed from 
the light curve. From top to bottom, the periods of the planet can- 
didates are 3.25 days, 6.20 days, 8.35 days, 13.77 days and 87.09 
days, respectively. The candidates with periods of 13.77 days and 
6.20 days have already been identified as KOI 505.01 and KOI 
505.02. 



date, really similar to KIC 009535585, with 0.19 day period 
and 1.5 R® radii. (See Figure[8|e), row 4, column 3.) 

KIC 005689351 (KOI 505) : A rich system with po- 
tentially as many as 3 more planet candidates. We initially 
identified a 2.2 ii® planet candidate with a period of 8.348 
(a) day in the system, which already has 2 KOIs. While this 
paper was under revision, 012 pointed out that there might 
be 2 additional signals in the system, with periods of 3.25 
days (b) and 87.09 (c) days. (See Figure [| (b), row 3, col- 
umn 1, 2, 3.) We find the same signals in this light curve 
after reprocessing the data. The folded light curve of the 
three planet candidates are separately shown in Figure [14] 
When modeling one planet candidate, we make use of the 
detected periods and epochs to filter out the transit signals 
due to the other two. The top and middle panels present the 
two KOIs, with periods of 13.767 days and 6.1956 days. The 
third planet candidate we found is in a 4:3 resonance with 
the latter. This is the first reported 3:4:6 resonance 3-planet 



candidate system. Rein et al. (20121 suggests that a reso- 



nance chain of three or more planets, in analog to Jupiter's 
moon system, would overcome the difficulties of forming the 
4:3 resonance in the traditional way. This system provides 
an interesting testing-ground for the theory. 

KIC 007595157 (KOI 246): A sub-earth candidate 
with 0.64 i?0 orbiting the star with a period of ~ 2.4 days. 
(See Figure[8](c), row 2, column 2.) It was a single planetary 
system in B12 catalog, with KOI 246.01 (1.3 R® radii and 
3.4 days period). 

KIC 008753657 (KOI 321): A candidate hot Mars, 
also the smallest candidate from our new findings, with a 
modelled radius of 0.5 R® and period of 4.6 days. (See Figure 
[8](c), row 6, column 1.) This is also reported by the Ofir team 



(012) with similar modelled parameters. The system was 
previously known to host KOI 321.01, a 2.4 days candidate 
superearth. 

KIC 11295426 (KOI 246): B12 catalog found it to 
host a superearth (KOI 246.01) with a 5.38 days period. We 
found an outer planet candidate with 0.58 i?® radii and 9.6 
day period. (See Figure[8](e), row 1, column 1.) This is also 
reported by 012. 



4 CONCLUSIONS 

We have analyzed 124840 stars with public Kepler data from 
quarters Q1-Q6 in total. The large majority of our candi- 
dates have already been identified by the Kepler team. We 
recover 92% of Kepler findings from the Bll catalog and 
the majority (86.4% of the planet hosting stars and 82.5% 
of the planet candidates) of the B12 catalog. Forty new false 
positives are identified in our analysis of centroids variation. 
We report 150 new planet candidates and 7 single transit 
events that haven't been assigned as KOIs (Kepler Objects 
of Interests) in this blind search. 55 of these planet candi- 
dates are listed as potential transit signals by the automatic 
TPS algorithm developed by the Kepler team (T12). While 
this paper was under revision, 012 conducted a search of 
new planet candidates using the SARS pipeline in all the 
KOIs, 43 of the planet candidates in this work overlap with 
their findings. To our best knowledge, 22 of these planet can- 
didates and 3 of the single transits are also independently 
identified by the Planet Hunters in their public website dis- 
cussion session. 40 of the planet candidates and 4 of the 
single transits are reported for the first time. 

The periods of our new candidates range from 
~ 0.17 days to > 440 days. The estimated planetary radii 
vary from ~0.5i?® to 62 R®. 124 of the planet candidates 
and three of the single transits have sizes smaller than 3 R® . 
We also found 6 new candidate multiple systems. In addi- 
tion, we found 57 more planet candidates in the already 
known KOI systems. By comparing our new findings with 
the Kepler candidates, it appears that for short periods 
(< 10 day), our pipeline and search procedure finds some- 
what more candidates than the Kepler pipeline. Most of the 
single transit events (and a few of the planet candidates) 
were found by visual inspection of the light curves that were 
flagged by BLS. This suggests that combining automated 
searches, and visual inspection is an efficient approach for 
transiting planet searches. Our searching process could po- 
tentially act as a supplementation of the findings from the 
Kepler team to improve the statistics in the short period 
planets. We trust that independent searches will be benefit 
both the Kepler team and the community. 
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Table 1. Anomaly Summary Table 



Start(BJD-2454000) 


End(BJD-2454000) 


Quarter 


Anomaly Type 


1002.5198 


1002.7241 


2 


EXCLUDE a 


1014.5146 


1016.7214 


2 


SAFE MODE b 


1033.2932 


1033.3341 


2 


TWEAK c 


1056.4853 


1056.8313 


2 


COARSE POINT d 


1063.2896 


1064.3726 


2 


EARTH POINT ° 


1073.3632 


1073.3632 


2 


ARGABRIGHTENING f 


1079.1866 


1079.1866 


2 


TWEAK 


1080.8008 


1080.8212 


2 


ARGABRIGHTENING 


1088.4018 


1089.3213 


2 


COARSE POINT 


1093.2239 


1093.2239 


3 


EARTH POINT 


1100.4163 


1100.5389 


3 


COARSE POINT 


11U4.0ZOO 


110/1 £/1QC 


Q 
O 




1113.4116 


1117.3346 


3 


MEMORY ERROR s 


1123.5461 


1124.4248 


3 


EARTH POINT 


1206.2584 


1206.2584 


4 


TWEAK 


1216.4137 


1217.3128 


4 


EARTH POINT 


1229.8386 


1233.8231 


4 


SAFE MODE 


1307.9995 


1309.2869 


5 




1336.7713 


1337.6091 


5 




1399.5042 


1400.3828 


6 


EARTH POINT 


1431.2374 


1432.1978 


6 


EARTH POINT 



a EXCLUDE: manually excluded cadence before pipeline processing. 

b SAFE MODE: due to unanticipated sensitivity to cosmic radiation, or unantici- 
pated responses to command sequences. The amplitude of flux is affected after the end 
date of safe mode, data for additional 3 days are removed rather than corrected. Data 
on both side are offset in magnitude to ensure continuity by a polynomial fit. 

c TWEAK (Attitude Tweaks): discontinuities in the data due to small attitude 
adjustments. The discontinuities are corrected by offseting the data on both sides using 
a polynomial fit. 

d COARSE POINT (Loss of Fine Pointing): due to losing fine pointing control, 
removed. 

e EARTH POINT: change of attitude due to monthly data downlink. Affect data 
the same way as safe mode. Corrected by the same method as safe mode. 

f ARGABRIGHTENING: diffuse illumination of the focal plane. Removed. Only 
argabrightening events longer than one cadence are listed here. 

s MEMORY ERROR: due to onboard spacecraft errors, gapped by the pipeline, the 
continuity on both ends of a gap is checked. Only a single memory error longer than 
one cadence is listed here. 

h The format of data release handbook 8 for Q5 data is different from other quarters. 
There are no safe mode, tweak, coarse point or exclude phenomena in Q5, only two big 
gaps, possibly due to earth point, are listed. 
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Table 3. False Positive List 









RA 








I CI lOLL 




uio^jideeiiieiiL y 


DIN J\ 














(R TD-94 c U00n^ 


\ ud J J 


(Pixel) 


(Pixel) 




2442359 


13.931 


19 


24 53 


.566 


+37 46 53.58 


965.8428 


0.55283 


-0.0025 




26.16 


3336765 


13.617 


19 


19 49 


.397 


+38 27 49.07 


965.1613 


1.84483 




-0.00032 


38.93 


3852258 


13.819 


19 


27 55 


.433 


+38 58 16.25 


968.1106 


5.75838 




0.005 


17.48 


4035667 


10.035 


18 


58 29 


.566 


+39 10 56.68 


965.5837 


2.87366 


-0.007 


0.005 


-1.00 


4072333 


15.664 


19 


41 19 


.236 


+39 09 45.14 


965.5638 


2.00775 




0.0030 


37.69 


4077901 


13.007 


19 


45 42 


.948 


+39 08 01.54 


967.7249 


6.05444 


+ 


0.002 


27.22 


4270565 


15.154 


19 


34 24 


.946 


+39 23 18.96 


972.2024 


18.01143 


0.007 


-0.007 


28.47 


5443775 


12.936 


19 


21 20 


.671 


+40 39 11.02 


967.8765 


3.30754 




-0.00034 


73.09 


5565497 c 


11.521 


19 


57 31 


.702 


+40 45 29.30 


966.11575 


1.412532 


-0.0126 


-0.0011 


131.4 


5622812 


12.832 


19 


30 41 


.270 


+40 52 10.34 


964.8384 


0.10304 


-0.005 


-0.005 


49.76 


5649325 c 


13.560 


19 


55 55 


.289 


+40 52 19.34 


1196.35 


195.63 


0.0016 




7.81 



a We only report the false positives not in the Kepler FP catalog here. We do not report FPs in single transit events. These FPs 
are all identified by examining the centroids. We refer to the text for a detailed description of the method. A complete version of 
this table can be accessed from the online version and downloaded from https:/ /sites. google.com/site/largepunkelephant/tables. 

b The flux centroid displacement amplitude is computed from the phase folded centroids. The x and y direction are listed separately. 
We do not list a displacement when the shift is not detected. 

c There are visible companion(s) in the 2MASS image stamp within 20" x 20". 
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Table 4. Stellar Parameter Table 



Star 




J eft 


l°g3 


fpp /HI 
[re/xlj 

(cgs) 




u a h 


u h c 


2985587 


13.910 


6023 


4.278 


-0.069 


1.276 





.34 


0. 


29 


3128552 


14.523 


5530 


4.673 


-0.248 


0.761 





.41 


0. 


25 


3240049 


11.557 


4435 


2.127 


-0.065 


17.188 





.65 


0. 


09 


3245969 


15.681 


4825 


4.789 


0.054 


0.594 





.59 


0. 


13 


3328026 


15.147 


5681 


4.508 


0.147 


0.947 





.42 


0. 


24 


3345675 


15.635 


4105 


4.628 


0.137 


0.598 





.57 


0. 


17 


3346154 


14.575 


5513 


4.490 


0.118 


0.957 





.45 


0. 


22 


3439096 


13.799 


4940 


3.119 


-0.340 


5.493 





.52 


0. 


18 


3541946 


13.597 


5537 


4.728 


-0.140 


0.711 





.43 


0. 


24 


3558849 


14.218 


5938 


4.432 


-0.410 


1.052 





.33 


0. 


30 


3728432 


15.646 


4371 


4.583 


-0.254 


0.696 





.59 


0. 


13 


3764879 


14.024 


5845 


4.641 


0.006 


0.809 





.38 


0. 


.27 


3834322 


15.397 


4627 


4.631 


-0.004 


0.698 





.62 


0. 


11 


4150804 d 


12.888 













.40 


0. 


25 



a A complete version of this table can be accessed 
from the online version of the paper and downloaded from 
https://sites.google.com/site/largepunkelephant/tables. 

b The quadratic limb darkening parameter a. 

c The quadratic limb darkening parameter b. 

d The limb darkening coefficient is obtained from the estimated stellar 
properties based on the J, H and K magnitudes. 
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